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AB-^TRACT 


An improved theory for the prediction of the rotordynamic coefficients 
of turbulent annular seals has been developed since the original, 15 February 
1983, report [l] on this project. This supplentary report compares predictions 
from the new theory to the experimental results of [1] and also introduces a 
new approach for the direct calculation of empirical turbulent coefficients 
from test data. 

An improved short-seal solution is shown to do a better job of calculating 
effective stiffness and damping coefficients than either the original short- 
seal solution or a finite-length solution. However, the original short-seal 
solution does a much better job of predicting equivalent added-mass coefficient. 
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INTRODUCTION 


In the original report [1] on this contract, experimental results 
were compared to a "short-seal" theoretical model [2] which employed 
Colebrook*s friction-factor formula [3] for predicting turbulent friction- 
factors. Since reference [1] was completed, an improved finite-length 
solution procedure [4]* has been developed, and the data have been reanalyzed 
to directly calculate the Yamada-HIrs [5,6] empirical coefficents. This 
supplementary report provides a comparison between the experimental data of 
[1] and the new theory of [4] with empirical friction factor-coefficients 
which have been directly obtained from the data. 


*The analysis is included as Appendix A. 
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ORIGINAL PAfj^ 

OF POOR QUA'IITY 


IDENTIFICATION OF EMPIRICAL TURBULENCE COEFFICIENTS FROM TEST DATA 

The finite-length-solution development [4] is provided in Appendix A, 
The leakage formula provided Eq. (15) of this reference is: 

^ p s /(ijLi) + [ 2 g - 2 cr 6 (Hmo) + 4q] ) p 

(0 (1 - ) 2 


Where 

p: Fluid density. 

^ = Entry loss coefficient. 

q = = Taper parameter. 

Co, Ci: Seal entrance and exit clearances, respectively, 
a = A L / C 

1+mo 

A = no R (1 + ) ; Wall friction factor, 

ao 4t)2 

Rao = 2 Y C / v: Centered-position Reynolds number. 

C" = (Co + Cl) / 2: Average seal clearance. 

L: Seal length. 

Y = Fluid Kinematic viscosity, 
b = V/RW 

Y = Q/ 27t R Average axial fluid velocity. 

Q = Volumetric flow rate. 

R = Seal radius, 
w = Shaft arfgular velocity. 

3 = 1 / (1 + 4b2) 

mo, no: Empirical coefficients for the friction-factor definition. 
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The data for each dynamic seal test includes the inlet and exit chamber 
pressures and five pressure measurements within the seal. Two of the pressure 
measurements within the seal are immediately interior to the inlet and exit. 
The volumetric flowrate and inlet and exit temperatures are also measured. 

Our objective is to take this data and det -X' rmine the entry-loss coefficient 
and the empirical friction"factor coefficients mo, no . 

For the ith test case, the total entry-loss factor (1 + Ci) is readily 
calculated from the inlet pressure drop relationship 


= ( 1 + gi ) • fiSL" 

(1 + q )2 2 


( 2 ) 


Eq. (2) was solved directly for (1 + ^i) for each test case, 

The calculations of mo, no for a given test case is relatively straight 
forward provided X can be determined from experimental data, since 


X (1 + = no [R (1 — ) ] 

4b^ 4b" 


•2 mo 


(3) 




In [X(l + ”) " ] = In (no) + mo In [R (1 + ] 

4b^ 4b^' 


Eq. (3) is linear in the parameters In (no) and mo. Hence, a least-square 
curve fit of all cases for a given housing-rotor combination will yield the 
desired data. 

The determination of X from a given data set, i.e., a given rotor-housing 
combination, represents the principal complication in executing this procedure. 
The pressure drop within the seal, AP|r, is defined as the difference between 
pressure measurements immediately interior to the inlet and exit of the seal, 
for case i, Eq. (1) can be expressed; 


(G. + 4q) pV. 

= _! - — 1 

(1 - 2 


(4) 
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Where 


« E 0 ^ [1 *• 3 0+mo) q^J 

The quantity may be readily determined from the experimental data, 
Hence, one can solve for as 


X^. « C G^7 2L' [1 - 3 (1+mo) q^J (5) 

where L' is the distance between the inlet and exit pressure tops. Direct 
solution of X^ from tnis equation is complicated, because mo on the right- 
hand side is unknown. This difficulty is resolved by an iterative procedure 
wherein an initial value for mo is guessed, which permits an initial calculation 
of the X^-'s for all cases in a data 5et. After the least-square solution yields 
an estimate for mo, rio, the procedure is repeated using the updated estimate, mo’. 
Convergence is rapid since the product 3 (1 + mo) q^ of Eq. (5) is generally 
much smaller than unity. 

Application of the above procedures yielded the results of Table 1. 


Case Housing Rotor no mo 


•! 

1 

1 

1 

0.20163 

-.2796 

2 

1 

2 

0.07106 

-.19691 

3 

1 

4 

0.00213 

.15089 

4 

1 

5 

0.00985 

.00980 

5 

2 

1 

0.07822 

-.1182 

6 

2 

4 

0.03831 

-.0162 

7 

2 

5 

0.04330 

-.0405 

8 

3 

1 

0.02300 

-.0377 

9 

3 

4 

0.00394 

.1448 

10 

3 

5 

0.00513 

.1118 


Table 1. Empirical coefficients for the forced seal data of reference [1] 
as determined from test data. 
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Of these ten cases, only tiui first two have the same d1 recti onany-homo- 
geneous surface-roughness treatment for both the rotor and housing. In the 
remaining cases, a circumferentially-grooved surface roughness treatment 
was inscribed in either the rotor, the stator, or both rotor and stator, 
Complete dimensions and surface-roughness measurements for the rotors and 
stators are provided in reference [1], For comparative purposes, Yamada's 
test results for smooth rotor and stators were no = 0,079; mo = -0,25. 


FORCE COEFFICIENT CALCULATION AND COMPARISON 


The "finite-length” solution procedure of reference [4] can be run in 
either a finite-length mode or in an improved short-seal mode. The data 
of Table 1 were used v/ith the (improved) short-seal and finite-length options 
of reference [4] to calculate radial and tangential force components for 
comparison to the tapered-seal test data of reference [1], Figures 1 through 
10 illustrate the results for the finite-length solution, while figures 11 
through 20 illustrate the results for the improved short-seal solutions. 

The data from these figures were used to calculate effective stiffness, 
damping, and acJded-mass coefficients. The results for the finite-length and 
improved short-seal solution options are provided in tables 2 and 3 respectively. 
The equivalent comparison between the original short seal theory and experimental 
results are provided in tables 4 and 5. 

A review of the results for aJJ. data sets shows no clear superiority for 
any procedure. As expected, the finite-length solution consistently predicts 
smaller values for the seal coefficients than either of the short-seal solutions. 

As earlier, the analyses only strictly apply to the first two data 
sets fr:r which the same directional ly-homogeneous surface roughness holds 
for both the stator and rotor. To compare the solution approaches for these 
data sets, the following least-square error calculations were made: 


EK = 

n 

E 

[1 • 

- ( KEX/ KTH 

] 

(Stiffness) 

i 

= 1 





EC = 

n 

l 

[1 ■ 

- (CEX/CTH)^ 

J 

(Damping) 

i 

= 1 





EM = 

n 

1 

Cl ■ 

• (MEX/MTH)^ 

] 

(Mass) 

i 

= 1 






where KEX and KTH are the measured and theoretical equivalent direct stiffness 
respectively, etc. 
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6. Fj^/A and Fg/A versus w for rotor 4, housing 2. Measured [1] and 
finite-length theoretical results [4]. 
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7. F^/A and F^/A versus oj for rotor 5, housing 2. Measured [1] and 
finite-length theoretical results [4]. 
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14. F^/A and Fq/A versus w for rotor 5, housing 1. Measured [1] and 
improved short-seal theoretical results [4]. 
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15. F^/A and F^/A versus w for rotor 1, housing 2. Measured [1] and 
improved short-seal theoretical results [4]. 
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16. F^/A and F^/A versus w for rotor 4, housing 2. Measured [1] and 
improved short-seal theoretical result [4], 
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17. Fj,/A and Fg/A versus oj for rotor 5, housing 2. 
improved short-seal theoretical results [4]. 
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and Fg/A versui o) for rotor 1, housing 3. Measured [1] and 
improved short-sea'' vheoretical res'.^lts [4]. 
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19. F^/A and Fg/A versus w for rotor 4, housing 3. Measured [1] and 
improved short-seal theoretical results [4]. 
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20. F^/A and Fg/A versus w for rotor 5, housing 3. Measured [1] and 
improved short-seal theoretical results [4]. 
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damping, and added-mass coefficients 


Of?IGfNAl PAGE 

Op POOR quality: 


TAPERED SEAL THEDRV SEAL 2 (TAPERED) 
UNITS ARE N/N 



KEF 


KEF 


KEF 

CEF 


CEF 


CEF 


(EXP) 


(THEORY) 

(EXP/THE) 

(EXP) 


(THEORY) 


(EXP/THE) 

HOUSING 1 

ROTOR 1 










RA=435, 000 

0. 8264E 

07 

0. 9511E 

07 

0. B689 

0. 2099E 

05 

0. 19I9E 

05 

1. 094 

RA»325, 000 

0. 3852E 

07 

0. 525SE 

07 

0. 7330 

0. 1624E 

05 

0. 1501E 

05 

1. 082 

RA^^ISO, 000 

0. 5826E 

06 

0 1241E 

07 

0. 4695 

2863. 


6561. 


0. 4364 

RA« 7D. 000 

-0. 2BS9E 

06 

0. 15B9E 

06 

-1. 818 

1369. 


4329 


0. 3162 

RA“ 30,000 

~0. 1894E 

06 

0. 1651E 

06 

-1. 147 

522. 1 


3617. 


0, 1443 

HOUSING l‘ 

ROTOR 2 










RA“ 30,000 

0. 2082E 

06 

-0. 2951E 

06 

-0. 7055 

1098. 


6076 


0. 1807 

RA= 75,000 

0. 3946E 

06 

-0. 561 7E 

05 

'■*7. 023 

1451. 


6050. 


0, 2398 

RA«1S0, 000 

0 U92E 

07 

0. 1133E 

07 

1. 052 

4122. 


6B65. 


0. 6004 

RA«335, 000 

0. 4745E 

07 

0, S640E 

07 

0. 8413 

0. 13B2E 

05 

0. 1609E 

05 

0. 85B9 

RAC450, 000 

0 S973E 

07 

0, 9450E 

07 

0. 9495 

0. 205BE 

05 

0. 2044E 

05 

1. 007 

HOUSING 1 

ROTOR 4 










RA= 50,000 

0. 5094E 

06 

0. 2701E 

06 

1. 886 

B37. 2 


1711. 


0. 4B93 

RA« 73, 000 

0, 7992E 

06 

0. 4103E 

06 

1. 948 

1647. 


2632. 


0 6258 

RA^ISO, 000 

0. 1900E 

07 

0. 1B60E 

07 

1. 022 

462B. 


5800. 


0. 7979 

RA«325, 000 

0. B0I3E 

07 

0. S9a6E 

07 

0. 8917 

0. 2081E 

05 

0. 1371E 

05 

1. 518 

RA»345, 000 

0. 9892E 

07 

0. 9419E 

07 

1. 050 

0. 2089E 

05 

0, 1472E 

OS 

1. 419 

HOUSING 1 

ROTOR 5 










RA* 50, 000 

-0. 4863E 

05 

0, 1695E 

06 

-0. 2869 

755. 3 


1B15. 


0. 4161 

RAi= 75, 000 

-0. 4094E 

05 

0. 3645E 

06 

-0, 1123 

15B5. 


2B21. 


0. 5619 

RA=150, 000 

0. e2S7E 

06 

0. 1519E 

07 

0. 5456 

4749. 


5910. 


0. 8036 

RAs32S, 000 

0. 7241E 

07 

0, 772BE 

07 

0. 9370 

0. 1627E 

05 

0. 1500E 

05 

1. 0B5 

RA«3B0, 000 

0. 9206E 

07 

0. 9935E 

07 

0, 9266 

0. 1935E 

05 

0. 1662E 

05 

1. 164 

HOUSING 2 

ROTOR 1 










RA* 75, 000 

0. 130BE 

06 

0, 5072E 

06 

0. 2579 

3530, 


3498. 


1. 009 

RA=125, 000 

0. 7659E 

06 

0. 1412E 

07 

0, 5424 

67S7. 


5752. 


1. ISO 

RA=150, 000 

0. 13S6E 

07 

0. 1952E 

07 

0. 7100 

9197. 


6879. 


1. 337 

RA--=275, 000 

0. 3912E 

07 

0. 6323E 

07 

0. 9350 

0. 21B7E 

05 

0. 1229E 

05 

1. 779 

RA=330, 000 

0. 0049E 

07 

0, BS45E 

07 

0. 9100 

0. 2343E 

05 

0. 1445E 

05 

1. 621 

HOUSING 2 

ROTOR 4 










RA*= 75, 000 

0. 1017E 

07 

0. 6449E 

06 

1. 577 

5586, 


3331. 


1. 677 

RA=150, 000 

0. 2039E 

07 

0, 261SE 

07 

0. 7979 

0. 1152E 

05 

7171. 


1, 606 

RA«270, 000 

0. 5061E 

07 

0. 7640E 

07 

0, 6624 

0. 2117E 

05 

0, 12E7E 

05 

1. 645 

HOUSING 2 

ROTOR 5 










RA= 75, 000 

-0, 5772E 

05 

0. 6502E 

06 

-0. BS77E-01 

2734. 


3896. 


0. 7017 

RA*125, 000 

0. 2427E 

06 

0. 1605E 

07 

0. 1512 

6638. 


6318. 


1 051 

RA=275, 000 

0. 3S42E 

07 

0. 7853E 

07 

0, 4892 

0. 2044E 

05 

0. 1354E 

05 

1. 510 

RA“290, 000 

0. 4905E 

07 

0. B373E 

07 

0, 5858 

0, 1851E 

05 

0. 1384E 

05 

1. 337 

HOUSING 3 

ROTOR 1 










RA=106, 000 

0. 6500E 

06 

0. 7782E 

06 

0. B353 

3732. 


4030. 


0. 9261 

RA=163, 000 

0. 1447E 

07 

0. 1925E 

07 

0. 7517 

6972. 


6152 


1. 133 

RA=303, 000 

0. 721 IE 

07 

0, 681BE 

07 

1. OSS 

0. 1827E 

05 

0. 1238E 

05 

1. 476 

RA=>355, 000 

0, 9S57E 

07 

0. 9160E 

07 

1, 076 

0, 2243E 

05 

0. 1461E 

05 

1. 535 

HOUSING 3 

ROTOR 4 










RA=106, 000 

0. 159SE 

07 

0. 1 190E 

07 

1, 343 

2746. 


4370. 


0. 6284 

RA=J.64, 000 

0. 2536E 

07 

0, 2B77E 

07 

0. BB15 

7136. 


6990. 


1. 021 

RA=2S8, 000 

0. 71B1E 

07 

0. 8448E 

07 

0, 8500 

0. 1942E 

05 

0. 1325E 

05 

1. 466 

HOUSING 3 

ROTOR 5 










RA=10o- COO 

0. 721 4E 

06 

0. 991 BE 

06 

0, 7274 

3819. 


4687. 


0. 8148 

RA-155, 000 

C. 1495E 

07 

0. 2142E 

07 

0. 6979 

7248. 


6888. 


1, 052 

RA=305, 000 

0 SI24E 

07 

0. S7S0E 

07 

0. 9253 

0, 1984E 

05 

0. 1394E 

05 

1, 423 


Table 4. Measure(d and original short-seal theoretical [2] predictions for 
effective direct stiffness and damping coefficients (Table 34, 
ref. [1]). 
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ORIGINAL PASS m 
OF POOR QUALin 


I 






r 


ft 




t 


TAPCRED BEAL 

THEORY 

SEAL 2 (TAPERED) 

UNITS ARE kg 




MEF 

MEF 

MEF 


(EXP) 

(THEORY) 

(EXP/THE) 

HOUSING 1 

ROTOR 1 



RAS435. 000 

3. 619 

4. 952 

0. 7309 

RA*=325/ 000 

19. 20 

7. 563 

2. 539 

RAalSO/ 000 

4. 77B 

6. 709 

0. 7122 

RA“ 75, 000 

7. 571 

10. 77 

0. 7027 

RA« 50,000 

6. 075 

10. 18 

0. 5965 

HOUSING i 

ROTOR 2 



RA» 50, 000 

2. 412 

25. 65 

0. 9402E-01 

RAW 75, 000 

3. 071 

17. 92 

0, 1714 

RA=>150, 000 

5. 012 

a. 489 

0. S903 

RAW335, 000 

9 460 

7. 808 

1. 211 

RA=450, 000 

-2 615 

6. 481 

-0, 4034 

HOUSING 1 

ROTOR 4 



RAW SO, 000 

-3. 935 

3. 522 

-1. 117 

RAW 75, 000 

-4. 75'’ 

4. 074 

-1, 168 

RAwl SO, 000 

2. 362 

2. 275 

1. 038 

RA“325, 000 

12. 22 

0. 4564 

26. 76 

RA=345, 000 

-7. 406 

3. 546 

- 2 . ose 

HOUSING 1 

ROTOR 5 



RA= 50, 000 

6 240 

4. 852 

1. 286 

RAW 75, 000 

9. 222 

4. 972 

1. 855 

RAwlSO, 000 

1 1. 84 

5. 480 

2, 160 

RA«325, 000 

-1. 466 

4. 107 

-0. 3569 

RA»380, 000 

-7. 613 

3. 668 

-2. 076 

HOUSING 2 

ROTOR 1 



RA= 75, 000 

20 86 

3, 540 

5. 892 

RA = 125, 000 

29. 9S 

2, 42S 

12. 35 

RAwl*30, 000 

24. 98 

-0. 5142 

-48. 58 

RAw275, 000 

24. 51 

1. 419 

17. 27 

RAw330, 000 

36. 44 

1, 333 

27, 34 

HOUSING 2 

ROTOR 4 



RAW 75, 000 

-3, 496 

1. 999 

-1, 749 

RAwlSO, 000 

19. 79 

-0. 7368 

-26, 86 

RAw270, 000 

-3. 745 

1 . 129 

-3. 316 

HOUSING 2 

ROTOR 5 



RA= 75, 000 

16. 21 

2 098 

7. 726 

RA=125, 000 

25. 60 

2. 913 

8. 7SB 

RA=275, 000 

7. 951 

0. 2855 

27. 85 

RAw290, 000 

-5, 850 

3. 145 

-1. 860 

HOUSING 3 

ROTOR 1 



RAwl06, 000 

10. S3 

3, 31.6 

3. 176 

RAwl63, 000 

19. 94 

3. 103 

6. 425 

RA=303, 000 

11. 81 

2. 147 

5. 499 

RAW355, 000 

4 242 

2. 908 

1. 459 

HOUSING 3 

ROTOR 4 



RA=106, 000 

-0. 9733 

1. 198 

-0. 8123 

RAW164, 000 

7. 962 

0. 5052 

15, 76 

RA=i:eo, 000 

9 343 

1, 526 

6. 122 

HOUSING 3 

ROTOR 5 



RA=106, 000 

13. 70 

2. 238 

6. 121 

RA=155, 000 

17, 72 

1, 270 

13. 95 

RAW305, 000 

10. 18 

1. 262 

8. 070 




Table 5. Measured and original short-seal theoretical [2] 
predictions for effective added-mass coefficients 
(Table 40, ref. [2]). 
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The results of these error calculations are presented in Table 6 below. 



EK 

EC 

EM 

Finite 

1.618 

2.14 

111.7 

Improved Short 

.042 

1.27 

78.8 

Original Short 

.049 

2.96 

4.45 


Table 6. Least square error calculations for the 
first two data sets of Tables 2 through 5. 


Obviously, trnnimum values of error are desirable. For prediction of effective 
stiffness and damping coefficients, the improved short-seal solution is seen 
to be the best. However, the original short-seal solutions is much better 
for calculating the equivalent added-mass coefficient, More specifically, 
measured added-mass coefficients are much larger than predicted by either 
the finite-length or the improved short-seal solution. 
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APPENDIX A 


FINITE-LENGTH SOLUTIONS FOR THE ROTORDYNAMIC COEFFICIENTS 
OF CONVERGENT-TAPERED ANNULAR SEALS 


D. W. Childs 

Mechanical Engineering Department 
Texas A&M University 
College Station, Texas 77843 


ABSTRACT 

A combined analytical -computational method is developed to calculate 
the pressure field and dynamic coefficients for convergent-tapered high-pressure 
annular seals which are typical of neck-ring and interstage seals employed in 
multistage centrifugal pumps. Completely developed turbulent flow is assumed 
in both the circumferential and axial directions, and is modeled by Hirs' 
bulk-flow turbulent-lubrication equations. Linear zeroth and first-order 
perturbation equations are developed for the momentum equations and continuity 
equation. The development of the circumferential velocity field is defined 
from the zeroth-order circumferential -momentum equation, and the nominal 
pressure-leakage relationship results from the zeroth-order axial -momentum 
equation. 

The first-order perturbation yields three partial differential equations 
which are reduced to three ordinary, complex, differential equations in the 
axial coordinate z. These linear equations are integrated to satisfy the 
boundary conditions, and define the pressure distribution due to seal motion. 
Integration of the pressure distribution defines the reaction force developed 
by the seal and the corresponding rotordynamic coefficients. The solution 
does not employ linearization with respect to the magnitude of the taper 
angle or the degree of swirl. Finite-length solutions are compared to 
"short-seal" solutions. 
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NOMENCLATURE 


b: 

c, C; 

^d* 

f(z): 

h(z) « H/C: 

hi = Hi/C; 


«W «N* 

k, K 
m, M 

mo = -0.25 
no = 0.079 



Pi: 

Pi* 

q: 


t: 


= U^/Rw 
Ug = Ug/RW 

'^ 00 ’ % 1 * 

*^Z0’ ^Zl* 
v: 


z = Z/L 


Dimensionless coefficients defined In Appendix C. 
Dimensionless coefficient defined In Eq.(9). 

Dimensionless damping coefficients defined by Eq.(34), 

Seal discharge coefficient defined by Eq.(15), 

Dimensionless clearance function defined by Eq.(2). 

Dimensionless clearance function. 

First-order perturbation clearance function defined 
by Eq.(4) and (18). 

Dimensionless seal stiffness coefficients defined by 
Eq.(34). 

Dimensionless mass coefficients defined by Eq.(34). 
Coefficients for Hirs' turbulent lubrication equations. 

Fluid pressure (F/L^). 

Zeroth^'order perturbation pressure Introduced In 
Eq.(4), (F/L^). 

First-order perturbation pressure Introduced in 
Eq.(4), (F/L^). 

Dimensionless perturbation pressure defined In Eq, (8). 
Taper-angle parameter defined In Eq.(3). 

Independent variable time (T). 

Dimensionless axial and circumferential velocity components. 

Zeroth and first-order perturbations In Ug. 

Zeroth and first-order perturbations in u^. 

Dimensionless swirl variable introduced in Eq.(7), 
and defined by Eq. (11) . 

Initial (z=0) swirl. 

Dimensionless axial coordinate. 



rr: 

Cq, C^: 
H(z,0,t): 

Hq(z): 


Nominal seal radial clearance, (L). 

Entrance and exit clearances, respectively, (L). 

Clearance function, introduced in Eq.(4), and defined 
in Eq.(17), (L). 

Centered-clearance function defined by Eq.(2), (L). 


Hi(0,t): 

L: 

P • 

AP: 

R; 

= p(Rw)H/p: 

R^ := 2pVH/y: 

Rco p(Rw)Cf/y: 

Rao " 2pVC7y: 

T - L/ V ; 

U = Rw: 

U^, U^: 

V(z): 

Y: 

X, Y: 

Z, R0: 

a: 

e: 

K 


First-order perturbation in H, introduced in Eq.{4), 

(L). 

Seal length (L). 

Seal supply pressure (F/L^). 

Nominal pressure-drop across seal (F/L^). 

Seal radius (L). 

Circumferential Reynolds number. 

Axial Reynolds number, 

Centered-position, circumferential Reynolds number. 

Centered-position, axial Reynolds number. 

Transit time for a fluid element to traverse the seal. 
Seal tangential velocity. 

Axial and tangential fluid velocity components (L/T), 

Centered-position axial fluid velocity (L/T). 
Centered-position average fluid velocity (L/T). 

Radial seal displacements (L). 

Spatial coordinates illustrated in Figure 2. 

Seal taper angle illustrated in Figure 2. 

Seal eccentricity ratio introduced in Eq.{4). 

Inlet pressure-loss coefficient. 


I 


Dimensionless friction-factor defined in Eq.(9). 

Dimensionless time. 

Shaft angular velocity (T"^). 

Shaft processional velocity (T”^), introduced in 
Eq.(23)_. 
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INTRODUCTION 


In a series of publications, Black et al. [1-3] have explained the consid- 
erable influence of seal forces on the rotordynamic behavior of pumps. Figure 
1 illustrates the tv/o seal types which have the potential for developing 
significant rotor forces. The neck or wear-ring seals are provided to reduce 
the back leakage flow along the front surface of the impeller face, while the 
interstage seal reduces the leakage from an impeller inlet back along the 
shaft to the backside of the preceding impeller. Pump seals are geometrically 
similar to plain journal bearings, but have clearance-to-radius ratios on 
the order of 0.005, as compared to 0.001 for bearings. Because of the clearances, 
and normal ly-experienced pressure differentials, fully-developed turbulent 
flow normally exists in pump seals. 

As related to rotordynamics, analysis of seals has the objective of defining 
the reaction forces acting on the rotor as a consequence of shaft motion. For 
small motion about a centered position, the relation between the reacuion-force 
components and shaft motion may be expressed by 


X 


K k 
-k K 


C c 
-c C 


M m 
-m M 


0 ) 


Unlike hydrodynamic bearings, seals develop significant direct stiffness in the 
centered, zero-eccentricity position due to the distribution of the axial pressure 
drop between (a) inlet losses and (b) an axial pressure gradient due to friction 
losses. Further, experiments [2] have shown that the above relationship holds 
for fairly large eccentricities on the order of 0.5; i.e., the dynamic coefficients 
(K,k,C,c,M,m) tend to be relatively insensitive to changes in static eccentricity 
ratios. 

Prior analysis to define seal rotordynamic coefficients has involved the 


following developments: 

(a) Black and Jenssen [2], [3] used a bulk-flow analysis, with the 
circumferential bulk-flow velocity assumed to be fully-developed shear flow 
at ^ . In these analyses, the axial ^momentum equation incorporates Yamada's 
[4] friction-factor results for flow through rotating concentric cylinders, 
with the friction factor defined by average circumferential and axial Reynolds 
numbers. In analogy to "short-bearing" solutions, a short-seal solution is 
developed, which accounts for the circumferential flow due to shear, but 
neglects that due to pressure. The short-seal solution provides a definition 
for the dynamic coefficients of Eq, (1). 

(b) In an appendix to [1], an approximate finite-length solution is developed, 
and correction factors are developed as a function of L/D ratios for the short- 
seal dynamic-coefficient solutions. 

(c) In [3], Black and Jenssen define the friction factor as a function 
of the local axial and radial Reynolds numbers, i.e., the local clearance. 

(d) Allaire et al . [5] used Black's model to numerically calculate 
dynamic coefficients at large eccentricity ratios. Further, while Black and 
Jenssen define seal coefficients in a coordinate frame that rotates at half 
the shaft angujar velocity, and employ a coordinate transformation to achieve 
stationary-reference results, Allaire et al . perform all calculations in a 
stationary reference frame. 

(e) Black et al . [6] combined prior seal -analysis governing equations 
with equations previously derived for the analysis of "Journal -bearings with 
high axial-flow in the turbulent regime," to examine the development of circum- 
ferential flow in a centered seal as a function of axial seal position. They 
demonstrate that the circumferential velocity starts from an arbitrary initial 
velocity and asymptotically approaches as it proceeds axially along the 


seal. Stated differently, they account for the influence of inlet sv/irl. 
Predictions of the stiffness cross-coupling coefficient are generally reduced 
if the development of circumferential flow is accounted for in seal analysis. 
This analysis does not include the dependence of the friction-factor on local 
Reynolds numbers, i.e., local clearance introduced in [3], 

(f) Childs [7] performed an analysis of straight turbulent seals for 
small motion about a centered position based on Hirs' turbulent lubrication 
equations [8], The short-seal analysis was employed under less restrictive 
assumptions than those previously employed to derive seal dynamic coefficients. 
A single derivation, from one set of governing equations, yields results which 
include all previous "short-seal" solution developments. 

(g) Childs [9] completed a finite-length solution for straight turbulent 
seals using the Hirs-ba$ed model of [7]. 

(h) Fleming [10] analyzed straight seals with one-step and convergent 
tapered seals; concluding that optimally tapered seals can develop considerably 
higher direct stiffnesses than straight seals. Fleming's analysis yields only 
the direct stiffness term, and does .not include the effect of swirl; hence, 

his results are not adequate for a rotordynamic analysis of pump response or 
stability. Childs [11] performed short-seal analysis of convergent-tapered seal 
based on Hirs lubrication equations which defines all of the required dynamic 
coefficients of Eq. (1). 

The present analysis yields finite-length solutions for convergent- 
tapered seal geometries. The model is analyzed using the method of reference 
[9]; however, unlike preceding analyses, linearization assumptions are not 
required with respect to the magnitude of either the taper angle or swirl. 
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Seal Geometry 

Figure 2 illustrates the seal geometry. The clearance at the centered 
position is defined by 

Hq = (IT + %) - aZ = ^■[1 + q(l - 2z)] = ff (2) 


where a is the seal taper angle, and 



c = 

(Cq 

+ C^)/2, Z 

= Z/L, q = ^ 

2C ^0 

Cl 


(3) 

The 

ratio 

of 

t'ntrance to 

exit clearances is 





'=0. 

1 + 

_9. 






Cl 

1 - 

q 





The 

clearance 

ratio Cn/C 
0 1 

is the following tabular 

function of q 


Cq/C 

1 

00 

7 

3 2 


T.67 

1.285 

q 


1 

0.75 

0.5 0.333 


0.25 

0.125 


where q = 1 corresponds to a zero-clearance exit. Given that Fleming's optimum 
stiffness choices for are between 1.8 and 2.2, maximum values for q to 

0 i 

be expected in practice would be less than 0.4. 
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Simplified Perturbation Equations 


Hits' governing equations are provided in Appendix A, and are thoroughly 
discussed in reference[8] . These bulk-flow equations define the axial and 
circumferential velocity components (u^jUq) and the pressure, p, as a function 
of the spatial variables (Rq,Z) and time, t. The equations are expanded in 
the perturbation variables 

u^ ” U^Q t £^^21 * H = Hq + eH-j^ (4) 

^6 " "^00 ’ P - Po + ePi 

Where e = e/C is the eccentricity ratio, to yield the perturbation equations 
of Appendix B. 

These perturbation equations may be markedly simplified by carrying out 


the following steps: 

(a) Introduce the following nondimensional variables 

z - Z/L, T - t/T (5) 

where T is the fluid transit time defined by 
T = L/V (6) 

(b) Introduce the swirl variable v defined by 

Uqo = ^2 + V (7) 

(c) Introduce the nondimensional perturbation pressure 

Pi = Pi/pV" (8) 

where V is the average fluid velocity. 

(d) Identify the friction-factor coefficient 

C = ^ = = ^ (9) 


The parameter X can be factored out of the terms noA^. occurring in 
Appendix B. This factoring step yields the a^. coefficients of the 
following equations, as defined in Appendix C. 



Following these steps, the governing equations become 
leAotk-OhdeA 'Equations 

(a) Axial “Momentum Equation 

i-.. 

(b) Circumferential “Momentum Equation 

4 . cr - - 


OF POOR QUALITY, 


( 10 ) 


( 11 ) 


EAA6t’-0^deA Eqimtlom 

The first-or(ier equations of AppentJix B are a(d(iitionally simplified by 
substitution from Eqs(lO) and (11) to yield 
(a) Axial “Momentum Equation 

9Pi hi r Oa.^ 9nl 

- ” “(l “ mo)cja — + oa + (1 + mo) — 5 — + — 

8 ^ ' 0.4 L 0 ' ' i. J 


UqI 

+(1 + mo)aa«--^ + 

^2f3 


1 

b . 


II 

bf^ 


( 1 ^) 


9U7i 9U7i 1 9U7i1 

+ Uil{h + v)-r^ + ^ 


3t 


80 


(b) Circumferential “Momentum Equation 

8 p, h 

- = “(1 “ mo)crai 

7bf 


TW 


“ = “(1 “ mo)craj-^ + [ 203 ^ + (1 + mo)aa^] - — ^ ( 13 ) 

r (1 + mo)aa 2 
L ?f2 2b^ 


•. “zi + f; 


1 

P“ei 

b 

8t 


+ ijiTih + v)- 


1 ®“ei] 


(c) Continuity Equation 


8 u 


Z1 


^(D^-^^zi = 4[^-t(..v) 


8 z 


8 hi 8 hii 


(14) 


In contrast to earlier developments [7, 9, 111, q and v are not treated as 


small parameters in obtaining these equations. 


Zeroth-Order Perturbation Solutions 


The zeroth-order continuity equation has the solution HqU^q *= constant, 
and the centered-position axial -velocity distribution is accordingly defined 
in terms of the volumetric flowrate Q and cross-sectional area by 

V(z) = Q/ZttRHq = Q/2iTRCf 7/f 
where V is the average or mid-seal velocity. Hence, 

Uzo = V(z)/Rw = b/f 

Eq.(ll), the circumferential -momentum equation which defines v, is non- 
linear, but may be integrated numerically without difficulty. Alternately, 
linearization of Eq.(ll) in terms of q and v yields a reasonable approximation 
of the nonlinear solution, [11], The nonlinear numerical solution is used in 
the present study. 

Linearization of the zeroth-order axial -momentum Eq.(lO) in q and v is 
helpful in providing an initial estimate for leakage, and yields the following 
steady-state relationship 

AH _ P pF pV^ I 1 + C , 4q + 2a[l - 3(1 + mo)_q^]_ . . 

AP - Cd -2 2- j JTT-^ + . (1 . q2)z (15) 

where 

3 = 1/(1 + 4 b ^), 


and 5 is the entry-loss coefficient. The term 


(1 + q)^ 


(16) 


accounts for the total pressure drop at the inlet, while the remaining terms 
account for the pressure drop due to wall friction and Bernoulli effects. For 
a specified AP, Eq.(15) may be solved iteratively for the average velocity V 
and associated leakage. The exact solution is obtained by iteratively solving 
the coupled differential Eqs.(lO) and (11). 




First-“Order Equation Solutions 

The preceding equations define Pj(z,0,t), and Ugj^{z,e,T) 

resulting from the seal clearance function The clearance H is 

defined in terms of the components of the seal-journal displacement vector (X,Y) by 
H = Hq - X cose - Y sine (17) 

Hence, by comparison to Eq. (4) 

ehj == -X cose - y sine (18) 

where 

X = X/C , y = Y/C 

Note that h^ is not a function of z, and its time dependency arises from the 
displacement variables x(t), y(t). 

Solutions for the equations cited above must satisfy the circumferential 
continuity conditions 

U 2 i(z>T,e) = u^^ (z,T,e+2Tr) 

Uq^(z,T, 6) = UQ^(z,T,e+2TT) 
p^(z,T,e) = p^(z,t,0+27t) 

To satisfy these conditions, the following solution format is assumed 
U 2 i(z>T,e) = U 2 -jq(z,t) cose + 

Ug^(z,T,e) = Ug^^(z,x) cose + Uq^ 5 (z,t) sine (19) 

Pj^(z,T,e) = p^^(z,t) cose + p-j 5 (z,t) sine 

Substituting from Eqs. (18) and (19) into Eq. (12) eliminates e as an independent 
variable, and yields two real equations. By introducing the complex variables 
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Appendix C; Definition of a| 

2a^B =6^+02 
a^B = (v + + (v - %)02 

“(f) C^'' * + (v - -j)G^] 

836 =63+64 

V “(f) ■*■ ^>^<53 + (v - 52)^64] 
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(c) Continuity Equation 
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Appendix B; Tapered Seal Pe rturbation Equations 

■ ml I III II ■ I I l | ^ lJ »_ i -» I I'll'l , 1 | | I, I f II | | . 

Substitution of the perturbation variables of Eq.(2) into the equations 
of Appendix A yields the following perturbation equations: 

Zeroth“Order Equations 

(a) Axial -Momentum Equation 

"^0^ ^ - no ol+nio A + D u u 
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(b) Circumferential -Momentum Equation 
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(c) Continuity Equation 


Uzo = b/f 


First-Order Equations 

(a) Axial -Momentum Equations 
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; H1rs* Turbulent Lubrication Equations 
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Hirs' turbulent lubrication equations [8] define a bulk-flow theory which 
does not explicitly make any assumptions concerning either (a) local flow 
velocity due to turbulence, or (b) the shape of average flow- velocity profiles. 
Only the bulk-flow relative to a surface or wall and the corresponding shear 
stress at that surface or wall are considered or correlated. Hirs' axial and 
circumferential momentum equations can be stated, respectively, as 


zHi iP. _ no D l‘**nio 

yU 


pP _ 1!^ D 
+ R/ 


1+mo 

0^(00^ + u^^) 2 + 


1+mo 


9u. 


j_l 

IT 0t" T" IF" Iz 


(A.l) 


-H^ 1 9p „ no o 1+mo 
yU R 30 " 2 T 


1+mo 
2 


1+mo 




+ R, 


3u 


U 3t R 30 3Z 


0 


(A. 2) 


with the bulk-flow continuity equation 
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where = Cr^ is the amplitude of seal motion. The components are expressed 
as a function of because for a given seal geometry' (L,R,C) and set of 
operating conditions {AP,w), the excitation frequency qT is the only independent 
variable. Stated-differently, Sq. (33) provides a frequency-response solution 
for the reaction force components. 

To calculate seal coefficients, a comparable statement of reaction-force 
components is developed from the following nondimensional statement of Eq. (1) 
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Substitution from Eq. (32) yields 


XF^(fiT) 

ttRAPRq 


K + c(nT) - M(f2T)^ - 


t: 0^^ 


(35) 


XFrt(f2T) ~ n 1 

irRAPR ^ ^ ~ C(CiT) - m(S2T)^ = f32(^)dz 

Hence, the dynamic seal coefficients (K,k,C,c,M,m) may be obtained by comparing 
the solution obtained by Eq. (33) with Eq. (35). More specifically, they are 
obtained by a least-square curve-fit of the solutions stated on the right-hand 
side of Eq. (35) 
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Dy namic Coefficient Definitions 

Having obtained the pressure-field solution of Eq. (30), solution for the 
dynamic coefficients will now be undertaken. The reaction-force components 
acting on the rotor due to shaft motion are defined by 

Fj^(t) = -eRL PjCOS0d0dz - -cRLpV^ PjCOS0d0dz 

Ky(t) = -eRL PjSin0d0dz - -cRLpV^ p^sin0d0dz 

From the last of Eq. (18), these integrals further reduce to 

F^(t) = -eRLupV^ q/ p^^dz ; Fy(t) = -cRLirpV^ 

The motion defined by Eq. (22) is orbital at the precessional frequency and 
radius R^. This statement may be confirmed by comparing the last of Eq. (19) 
with Eq. (22) to obtain 


X - Cr^ cosQt , Y = Cr^ sin^t (32) 

Definition of the reaction forces is simplified by performing the integration 
at a time when the rotating displacement vector is pointing along the X axis, 
i.e., when fit = 0. Eq. (24) shows that and coincide for this time and 
location. Hence, Eq. (31) yields the following component force definitions 
parallel and normal to the displacement vector 


F^(JJT) = - r^(TrRLpV = ) f3c.(z)dz 

Fg(!^T) = - r^(7tRLpV2) y f3s(z)dz 
A useful nondimensional version of these equations is 
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A similar operation on Eqs.(13) and (14) yields 
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The time-dependency in these equations is eliminated by assuming a harmonic 
seal motion of the form 
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h^ = ~ = r^e'i^'^”^ 


(23) 


where r^ is a real constant. The associated harmonic solution can then be 
stated 
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P^(z,t) = p^(z) e 
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Substitution from Eqs. (23) and (24) into Eqs. (19) and (20) yields 
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where 
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-(1 “ mo)aa^/2f^ 

-[2q + (1 - mo)ga^]/f‘* + jTT/f^ 


The following three boundary conditions are specified for the solution of 
Eq. (25): 

(a) The exit perturbation pressure is zero, i.e., 

P^a) = 0 (27) 

(b) The entrance circumferential velocity perturbation is zero, i.e., 

Uq^(O) = 0 (28) 

(c) The pressure loss at the seal entrance is defined by 

- p(0,e,x) = I (O,0 ,t)(1 + 0 

This equation yields the following perturbation-variable boundary condition 


Pj (0) = - (1 + f.) + q) 
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Solution of the differential Eqs. (25) in terms of the boundary conditions 
Is relatively straightforward, yielding a solution for the velocity and pressure 
fields of the form 
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